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We point out the opportunity of tau neutrino astronomy for neutrino energies of the order of 
10 GeV to 10 3 GeV. In this energy range, it is demonstrated that the flavor dependence in 
the background atmospheric neutrino flux leads to drastically different prospects between the 
observation of astrophysical muon neutrinos and that of astrophysical tau neutrinos. Taking 
the galactic-plane neutrino flux as a targeted astrophysical source, we found that the galactic- 
plane tau neutrino flux dominates over the atmospheric tau neutrino flux for neutrino energies 
beyond 10 GeV. Hence the galactic-plane can in principle be seen through tau neutrinos with 
energies greater than 10 GeV. In a sharp contrast, the galactic-plane muon neutrino flux is 
overwhelmed by its atmospheric background until the energy of 10 6 GeV. 



1 Introduction 

The atmospheric v„ — > u T oscillations established by the Super-Kamiokande (SK) detector en- 
sures that a non-negligible u T flux reaches the Earth. The updated SK analysis of the atmospheric 
neutrino data gives^ 

1.9 • 1(T 3 eV 2 < 8m 2 < 3.0 • 1(T 3 eV 2 , sin 2 29 > 0.9. (1) 

This is a 90% C.L. range with the best fit values given by sin 2 29 = 1 and 5m? = 2.4- 10~ 3 eV 2 re- 
spectively. The tau neutrinos arising from the above — > v T oscillations are presently identified 
on the statistical basis^. On the other hand^, the total number of observed non-tau neutrinos 
from various detectors are already greater than 10 4 with energies ranging from 10" 1 GeV to 10 3 
GeV. It is essential to develop efficient techniques for the identification of tau neutrinos. 

For E > 10 6 GeV, the v T can be detected by the double bang show ers o f the tau lepton, which 
is produced by the neutrino nucleon charge current (CC) scatteringEEl In this case, the two 
showers are separated by the tau lepton decay length, which is roughly 50 m at E T = 10 6 GeV. 
For E < 10 6 GeV, one relies on the distinctive properties of the tau neutrino induced showers 
to detect v T . At low energy, the showers produced by CC v T interaction and the subsequent 
tau decay practically coincide. Treating these two showers as one, it is found that a much 
higher fraction of tau neutrino energy deposits in the form of showers than either in the 
CC interactions or in the neutral current (NC) interactions of any neutrino flavor. Hence, the 
footprints of v T might be identified by the study of energy spectrum of shower events^. The 

ratio of the combined electromagnetic and hadronic shower event rate to the muon event rate 
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is also a sensitive probe to tau neutrino appearance due to oscillations ' . The identification of 
v T appearance using shower properties requires that the flux of v T is comparable to those of 



other neutrino flavors. It is a challenging task to identify v T if other neutrino flavors dominate 
in flux. Nevertheless, the tau neutrino astronomy will be rather promising if the above tasks 
can be carried out, as we shall demonstrate in the following. 



2 The Prospects of Tau Neutrino Astronomy 

We illustrate the opportunity of tau neutrino astronomy by the possible detection of astrophys- 
ical v T from the galactic plane direction. To do this, we estimate the flux of galactic u T and 
calculate the background atmospheric v T flux in detail. The tau neutrino flux arriving at the 
detector on Earth, after traversing a distance L can be written as 

<fi?(E) = P{E) ■ <P^(E) + (1 - P(E)) ■ <p VT (E), (2) 

where (p Uti (E) and 4> Vt {E) are intrinsic muon neutrino and tau neutrino fluxes respectively, 
P{E) = Piyn — > v T ) = sin 2 28 • sin 2 (L/L osc ) is the ia, — > v T oscillation probability with the 
oscillation length given by L osc = AE/5m 2 . 



2.1 The galactic-plane tau neutrino flux 

One calculates the intrinsic galactic-plane Vn and v T fluxes by considering the collisions of inci- 
dent cosmic-ray protons with the interstellar medium. The fluxes are given by 

ME) = Rn p J™dE p <p p (E P ) da Y E " +¥ , (3) 

where E p is the energy of incident cosmic-ray proton, do~ pp ^, u+ Y/dE is the neutrino energy 
spectrum in the pp collisions, R the typical distance in the galaxy along the galactic plane, 
which we take as 10 kpc (1 pc ~ 3 • 10 16 m). The density of the interstellar medium n p along 
the galactic plane is taken to be ~ 1 proton per cm 3 . The primary cosmic-ray proton flux, 
(t)p(Ep) = dN p /dE p for E p < 10 4 GeV, is given by^l 

<t> p {E p ) = 1.49 • (E p + 2.15 • exp(-0.21^£ p )J , (4) 

in units of cm -2 s -1 sr -1 GeV -1 . 

The flux of galactic-plane muon neutrino arises from the two-body 7r decays and the sub- 
sequent three-body muon decays. The differential cross section for p + p ^ tt + X is model 
dependent. We adopt the parameterization in^for such a cross section. 

The galactic-plane v T flux arises from the productions and decays of t he D * mesons. It has 
been found to be rather suppressed compared to the corresponding i/„ fiux^. Clearly the total 
galactic-plane tau neutrino flux, (f>^{E), is dominated by the — > v T oscillations indicated by 
the term P(E) ■ <f> v JE) in Eq. ©. We have^ 



r(^) = 2-10- & ^J , (5) 
in units of cm~ 2 s _1 sr _1 for 1 GeV < E < 10 3 GeV. 

2.2 The atmospheric tau neutrino flux 

We follow a semi-analytic approach^ for computing the flux of intrinsic atmospheric Ua which 
could oscillate into v T . For 7r-decay contribution, the flux in the notation off*^ reads: 

d 2 NZ(E,e,X) r°° f E » ®(En-r^r) f x dX' 

1 dE N d E n -^—^4 ¥-P v (E v ,X,X') 



_ 5 f E yu* 



dEdX Je Je d n E„(l - j n ) Jo Xn 

( X'\ 



x—F Nw (E w ,E N )exp(-—)<l> N (E N ). (6) 
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Figure 1: The comparison of galactic- plane and the atmospheric v T fluxes in the presence of neutrino oscillations. 
The atmospheric v T fluxes in the left figure are the downward, while those in the right figure are the upward. 



We only consider the proton component of 0jv, which is given by Eq. @. The kaon contribution 
to the atmospheric v„ flux can be computed in the similar way. We note that charmed-hadron 
decays also contribute to the atmospheric v„ flux. However, these contributions are negligible 
for E < 10 5 GeV. 

The intrinsic atmospheric v T flux is reliably calculable using the perturbative QCD 12 . One 
writes the flux as 

d 2 N u ^(E,X) _ Z pDa Z DsUT exp(-X/A p )cf> p (E) 



dEdX 1 - Zpp(E) A p 

where the Z moments on the RHS of the equation are defined by 

(Ei ) Xi (Ej ) dn iA ^ jY (Ei , Ej ] 



(7) 



with driiA^jY(Ei, Ej) = daiA^jy(Ei, Ej) / OiA{E; L ). We have calculated the total atmospheric 
v T flux by applying Eq. © with </> v (E) given by d 2 N l/fir (E,X)/dEdX and integrating over 
the slant depth X. For £ < 70°, the oscillation probability P{yp, — > v T ) is calculated using the 
relation X = Voexp(— Lcos£//io)/ cos£, with Xq = 1030 g/cm 2 , ho = 6.4 km, and L the linear 
distance from the neutrino production point to the detector on Earth^l 

The comparison of the galactic-plane and the atmospheric v T fluxes is given in Fig. ^ The 
plot on the left hand side compares the galactic-plane v T flux and the downward atmospheric 
background flux. For 5m 2 = 2.4 • 10~ 3 eV 2 , sin 2 26 = 1, we find that both fluxes cross at 
E = 2.3 GeV. It is seen that the atmospheric v T flux is sensitive to the value of Sm 2 for E < 20 
GeV. This flux also changes its slope at E ~ 20 GeV. Below 20 GeV, the atmospheric v T 
flux predominantly comes from the oscillations, i.e., (/^(E) ~ 4> v {E) ■ sin 2 29 ■ s'm 2 (L/L osc ) 
following Eq. ©. Since L osc = 4^/<5m 2 w 330 km for E = 1 GeV and 5m 2 = 2.4 • 10" 3 eV 2 , we 
approximate sin 2 (L/L osc ) with (L/L osc ) 2 , so that ^"^(E) ~ 4> U)1 (E)E~ 2 . Because the neutrino 
oscillation effect steepens the § Vt spectrum for E < 20 GeV, the slope change of 4> Vt at E f» 20 
GeV is therefore significant. The plot on the right hand side of Fig. ^ compares the galactic- 
plane v T flux and upward atmospheric background flux. The two fluxes cross at E = 8 • 10 2 GeV 
for the best-fit neutrino oscillation parameters. We have also compared the galactic-plane and 



atmospheric v T fluxes for several other zenith angles. For instance, for the zenith angle £ = 60°, 
these two fluxes cross at E = 6.0 GeV for 5m 2 = 2.4 • 10~ 3 eV 2 with the maximal mixing. 

The above comparison of galactic-plane and atmospheric v T fluxes indicates a window of 
opportunity for the tau neutrino astronomy. Clearly, if v T can be identified from the overwhelm- 
ing z/„ background (in this case the atmospheric i>u flux), the galactic-plane can be seen through 
GeV energy tau neutrinos in the downward directions. In the upward directions, galactic-plane 
tau neutrinos are observable for E > 10 3 GeV. 

It is instructive to also compare the galactic-plane and the atmospheric Vu fluxes. While 
the galactic-plane neutrino flux is flavor independent, the atmospheric i/„ flux dominates its v T 
counterpart. As a consequence, the crossing point of galactic-plane a nd atm ospheric fluxes 
is pushed up to 5 • 10 5 GeV, which is drastically different from u T case ^^^ l This is a general 
situation in the astronomy where the opportunity for observing astrophysical neutrinos begins 
typically at 10 6 GeV. In contrast, with the future development of v T identification techniques, 
the energy threshold for the neutrino astronomy might significantly be lowered down. 
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